An accelerator mode based technique for studying quantum chaos 
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We experimentally demonstrate a method for selecting small regions of phase space for kicked rotor quantum 
chaos experiments with cold atoms. Our technique uses quantum accelerator modes to selectively accelerate 
atomic wavepackets with localized spatial and momentum distributions. The potential used to create the accel- 
erator mode and subsequently realize the kicked rotor system is formed by a set of off-resonant standing wave 
light pulses. We also propose a method for testing whether a selected region of phase space exhibits chaotic or 
regular behavior using a Ramsey type separated field experiment. 
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Chaos in quantum mechanics is still a relatively poorly de- 
fined concept. Typically it is taken to refer to the behavior of a 
quantum system which in the classical limit exhibits an expo- 
nential sensitivity to initial conditions. A much studied exam- 
ple of such a system is the delta-kicked rotor [JxL 0]. A realiza- 
tion of its quantum version, in the guise of cold atoms exposed 
to the ac-Stark shift potential of an off -resonant standing wave 
of light [ph, has elucidated many of the concepts associated 
with quantum chaos [[|, ||, ^, ^J. However, one substantial 
problem remains: there is no clear way of distinguishing reg- 
ular from chaotic dynamics. In the quantum regime it is not 
possible to define a chaotic region of phase space as being one 
in which two initially similar states have an overlap which de- 
creases exponentially with time. The unitary nature of any 
interaction necessarily implies that the overlap between such 
states remains unchanged. It has been suggested by Peres 
|§] that an alternative definition is needed for quantum sys- 
tems. Peres' proposal involves examining the evolution of a 
state which can interact with two potentials which differ very 
slightly in form. If the overlap between states produced by 
an evolution under each potential decreases exponentially as 
a function of time, then the region of phase space occupied 
by the initial state can be said to be chaotic in the quantum 
sense. The justification for this definition is that perturbation 
theory only converges when the region under examination has 
regular properties [^|. More recently Gardiner et al. JToj ] have 
proposed an experiment with a trapped ion which would im- 
plement the essential features of Peres' idea. 

In this paper we address two problems related to the exper- 
imental determination of a phase space stability map for the 
atom optical version of the QDKR. Firstly, we demonstrate a 
method for preparing cesium atoms in a restricted region of 
phase space using the same standing wave light pulses which 
create the rotor. Secondly, we show how by examining the 
overlap between the wavefunctions of the two hyperfine lev- 
els of the 6 2 Si/2 cesium ground state it should be possible to 
test the type of dynamics exhibited by the prepared atoms. 

The basis of our atom optical version of the QDKR is expo- 
sure of laser cooled cesium atoms to pulses (duration t p and 
separation time T) of a standing wave of off-resonant light. 
The pulses are short enough to allow the effect of the atoms's 



kinetic energy to be neglected during a pulse. This places our 
experiment in the Raman-Nath regime in which the spatially 
periodic ac-Stark shift potential created by the light acts as a 
thin phase grating for the atoms [jilt]. Thus an incident plane 
de Broglie wave is diffracted into a series of "orders" sepa- 
rated in momentum by TiG, where G = 2k\, (^l is the light 
wavevector). One of the striking features of this system is the 
existence of quantum resonances [f§|l§]. These resonances 
occur when the pulse interval is a multiple of the half-Talbot 
time, Ti = 2irm/frG 2 . During these special times all of the 
diffraction orders formed from incident plane waves with cer- 
tain momenta will freely evolve phases which are multiples of 
2ir. For cesium, the first quantum resonance occurs at T = 67 
//s. 



In addition to being a paradigm of experimental quan- 
tum chaos, this system can also produce quantum accelera- 
tor modes [|llj. This is achieved by adding a potential of 
the form Ua(x) = max, where m is the atomic mass, a is 
an applied acceleration and x is the position along the stand- 
ing wave. In our experiment the standing wave is oriented in 
the vertical direction, so a is the acceleration due to gravity. 
Quantum accelerator modes are characterized by a fixed mo- 
mentum transfer of (on average) £fi,G during each standing 
wave pulse (see Fig. 1 of Ref. [pi]]). Typically the accelerator 
modes are formed when the value of T is near to a quantum 
resonance. For simplicity, we will henceforth confine our at- 
tention to pulse repetition times near to the first quantum res- 
onance. One way of modelling the accelerator modes is to 
make the approximation that they consist of just a few diffrac- 
tion orders (say q — 1, q and q + 1, where q is an integer) which 
in the time between two light pulses accumulate a phase dif- 
ference which is very close to an integer multiple of 2tt. At 
the next pulse this makes it possible for interference to occur 
between the diffraction orders in such a way that the popula- 
tion of the three orders centered on q + £ [JHI] is enhanced. By 



using this rephasing condition, it can be shown [ 14 1 that after 



N p pulses the central momentum of the accelerator mode is 



given by q = -^-j^ — -^p£> where 7 = h^G 3 /2irm 2 a and 
a = T/Ti. The same rephasing condition also leads to the 
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conclusion that only incident plane waves with momenta 

Pinit = + const) TiG 



(1) 



(a) 



(where I is any integer) can ever participate in an accelerator 
mode. Thus in momentum space an accelerator mode resem- 
bles a comb with a tooth spacing of a~ 1 hG. Furthermore, 
since plane waves separated by this momentum behave iden- 
tically under the action of the kicks, the accelerator mode 
effectively contains only a single momentum. Figure [j](a) 
shows the theoretical momentum distribution of the acceler- 
ator mode after one set of pulses, as calculated using a model 
based on diffraction jilt]. To reflect the experimental situa- 
tion the starting distribution was gaussian with a width of 12 
hk^ at FWHM. Although there is good qualitative agreement 
between Eq. ([!]) and the numerically derived momentum dis- 
tribution, an explanation for the finite widths of the comb el- 
ements is needed. To understand this effect recall that each 
comb element must give rise to a set of diffraction orders 
which are always spaced by exact multiples of fiG. To de- 
termine what the real momentum distribution looks like we 
must add together the diffraction from all the different comb 
elements. These have a spacing of a~ 1 hG, which for pulse 
interval times just less than the Talbot time is slightly greater 
than TiG. The most obvious effect of including the diffraction 
orders is that the width of each comb element of the resultant 
distribution becomes non-zero and increases with a — 1 1 . Ad- 
ditionally, when all the diffraction orders are weighted by the 
intensity of the comb element from which they originated, the 
nearest-neighbor spacing of the peaks in the momentum dis- 
tribution is reduced as one moves away from the center of the 
distribution. The opposite effect is produced when the pulse 
interval is slightly greater than the Talbot time. Although the 
change in spacing is not immediately obvious in Fig. [jja), it 
has been confirmed in a detailed analysis of this data. 

In Fig.[j](b) we show the spatial form of the accelerator 
mode wavefunction as calculated with our numerical model. 
Such a distribution is similar to that deduced from the assump- 
tion that the mode consists of the sum of three plane waves. 
Importantly, the wavefunction is periodically localized in po- 
sition with maxima occurring every A spa t = 2tt/G. Since 
this is the wavelength of the standing wave potential, points 
having this separation behave equivalently and the wavefunc- 
tion effectively has a spread in position which is less than 
Aspat- From Fig.[j] we expect the width of one element of 
the momentum comb to be approximately 0.4 TiG and the ex- 
tent of each region of strong spatial localization to be A spa t/3. 
For comparison, the extent in momentum of a unit cell of the 
classical phase space is (Ti/ 2 /T)hG while that in position is 
Aspat- Hence the accelerator mode isolates a restricted region 
in phase space which is about 10% of the overall area. Addi- 
tionally, since the effect of the accelerator mode is to produce 
a large momentum offset, it should be straightforward to iso- 
late these atoms. 

We now discuss the experimental observation of the dis- 
tributions described above. Our atomic source was a cesium 
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FIG. 1 : Numerically simulated momentum (a) and spatial (b) distri- 
butions of the accelerated component of an ensemble of cold cesium 
atoms after 20 pulses of an off-resonant standing light wave. The 
pulses had a separation of 60 /is. The momentum distribution con- 
sists of a comb of peaks separated by a little more than two photon 
recoils. The degree of spatial localization is consistent with an accel- 
erator mode which contains approximately three diffraction orders. 
The dashed line in (b) shows the position of the standing light wave. 



magneto-optic trap containing approximately 10 7 atoms. Af- 
ter the trap was switched off the atoms were cooled in an op- 
tical molasses to a temperature of 5 /iK, corresponding to a 
momentum width of 12 Tik^ at FWHM. Following release, 
the atoms were exposed to a series of pulses from a verti- 
cally oriented standing light wave, detuned 30 GHz below 
the F = 4 — ► F' = 3 line in the Dl transition at 894 nm. 
The pulses had a duration of t p = 0.5 [is and a peak inten- 
sity of approximately 20 W/cm 2 . The atoms then fell freely 
until they passed through a probe laser beam resonant with 
the D2 cycling transition. By measuring the absorption of the 
probe light as a function of time we were able to determine 
the momentum distribution of the atoms in the F = 4 level. 
Additional details of our apparatus can be found in Ref. [0. 

The momentum resolution of our experiment was approxi- 
mately 2 photon recoils, determined by the initial spatial ex- 
tent of the optical molasses and the thickness of the probe 
laser beam. Hence direct observation of the accelerator mode 
momentum comb was not possible. Instead we used two sets 
of temporally separated light pulses to infer its existence. The 
first set contained 20 pulses and was used to create an ac- 
celerator mode. The resultant atomic distribution was then 
translated in momentum by allowing the atoms to fall for a 
variable amount of time t wa it ■ Finally, a second set of pulses, 
identical to the first, was applied. Figure ^| shows the resultant 
momentum distributions as a function of i wa it- Both panels 
contain a large fraction of unaccelerated atoms near p = 
(dotted line), a group of atoms that have been accelerated by 
one set of pulses (dot-dashed line), and atoms that have been 
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accelerated by both pulse sets (dashed line). In Fig. ^a) the 
pulses had an interval of 60 /is and the accelerator mode im- 
parted momentum in the negative direction (with the conven- 
tion gravity is negative), while in Fig.^(b) the interval was 
74 //s and the accelerator mode imparted positive momentum. 
The effect of t wa j t is to allow gravity to translate the distri- 
bution in momentum space. If any of this distribution over- 
laps with the momentum comb of Eq. ([[]) when the second 
set of pulses occurs then a further acceleration takes place. 
The periodic variation in the doubly accelerated population 
with t wa j t is just the length of time required for gravity to ac- 
celerate the atoms by the momentum separation of adjacent 
comb elements, a~ 1 TiG. In the T=60 /is case this gives a the- 
oretical value of 753 /is between accelerator mode revivals. 
Along the dashed line of Fig.^a) we observe 763 ± 4 /is, the 
discrepancy with the calculated value most likely due to the 
slight variation in comb spacing across the accelerator mode 
discussed previously. A similar level of agreement is found 
for T=74 /is. As regards the width of each momentum comb 
element, we note that numerical simulations have provided 
almost identical results to those shown in Fig.|| so it is rea- 
sonable to assume that the actual extent of each momentum 
element is as shown in Fig. |l](a). 

Although we have accounted for the accelerator mode re- 
vivals, we have said nothing about their sloping orientation 
(as seen in Fig. ||), nor about this slope's dependence on pulse 
interval. An explanation for this effect can be found by re- 
turning to the earlier observation that the momentum inter- 
vals between the peaks of Fig. |l](a) change as one moves away 
from the center of the accelerator mode. This implies that not 
all parts of the distribution have to be translated by the same 
amount to regain the accelerator mode condition of Eq. (jl|). At 
60 /is, the outer peaks bulge towards the center of the distribu- 
tion. Assuming that the whole distribution is accelerating un- 
der gravity, the first part to regain the accelerator mode condi- 
tion will be the component with the greatest positive momen- 
tum. For 74 /is the opposite will be the case and the direction 
of the revival's slope will flip. 

When £ wa j t is scanned in much smaller time steps the pre- 
vious experiment can also be used to demonstrate localization 
in position. These experiments are only meaningful if they are 
performed when the wavefunction has some degree of spatial 
localization, that is when i wa ; t is close or equal to an integer 
multiple of the pulse separation time. Furthermore, for these 
wait times all the elements of the momentum comb described 
by Eq. ([]]) translate by an integer number of standing wave 
periods (or equivalently acquire a phase of 2ir x integer). As 
we have seen, the momentum distribution produced by the ac- 
celerator mode is not a comb of delta functions. Thus we do 
not expect an exact rephasing of the distribution even at the 
special wait times. Experimentally we have found that for 
wait times greater than three pulse separation times very little 
rephasing is observed. Figure || shows the momentum dis- 
tributions observed when the pulse separation time is 60 /is 
and the wait time is scanned near to 120 /is. The accelerator 
mode reappears at the dashed line whenever the atoms pre- 




FIG. 2: Experimental momentum distributions of an ensemble of 
cold cesium atoms after two sets of 20 pulses from an off-resonant 
standing light wave. There was a variable amount of time t wa it be- 
tween the pulse sets. The dotted, dot-dashed and dashed lines denote 
atoms accelerated by zero, one and two sets of pulses respectively. 
The pulse spacing was T — 60 /is in (a) and T = 74 /is in (b). Each 
vertical slice is the momentum distribution obtained for one value of 
twait where the degree of shading indicates the population. Gravity 
is in the negative direction. The insets plot the data along the dashed 
lines (atoms accelerated by both pulses). 

pared by the first set of pulses have moved a whole number of 
standing- wave wavelengths. We can quantify what the period 
of the reappearance or revival should be by determining the 
distance which the atoms move during i wa it- This distance is 
given by Ax = vot wait + ^ai^ ait , where v is the velocity 
of the atoms after the end of the first set of pulses. We have 
calculated this velocity to be 177 mm/s fll5j], giving a revival 
period of 2.53 /is. This is in good agreement with the value 
of 2.56 ± 0.03 /is that is observed along the dashed line in 
Fig.||. Another important feature of this figure is that except 
for the exact multiple of the pulse separation time, not all parts 
of the accelerator mode produced by the first set of pulses are 
able to simultaneously undergo acceleration by the second set 
of pulses. This can be explained by noting that the distance 
moved by an atom during the wait time depends upon its ini- 
tial velocity. Thus the revival period for atoms that have been 
accelerated the most is the smallest. This picture also pro- 
vides an alternative way of understanding why we only see 
clear signatures of spatial localization near the multiples of T; 
the revivals from the different comb elements in the accelera- 
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FIG. 3: Experimental momentum distributions of an ensemble of 
cold cesium atoms after two sets of 20 pulses from an off-resonant 
standing light wave. Within each set the pulses were separated by 
T — 60 fis and there was a variable amount of time fwait between 
the sets. The dotted, dot-dashed and dashed lines denote atoms accel- 
erated by zero, one and two sets of pulses respectively. Each vertical 
slice is the momentum distribution obtained for one value of t wa it 
where the degree of shading indicates the population. The inset plots 
the data along the dashed line (atoms accelerated by both pulses). 



tor mode all come into phase at these times. 

The preceding results demonstrate that by using a pulsed 
periodic potential it is possible to select restricted regions of 
phase space. We now suggest a technique for investigating 
the dynamical properties of these regions under the action of 
a chaotic potential. Our proposal is to take two atomic states 
which have been prepared by an accelerator mode to have 
identical external components. We would expose each state 
to chaotic potentials which differ slightly in strength and then 
measure the overlap between the states in order to determine 
the dynamics. The F — 3 and F = 4 ground state hyperfine 
levels of cesium offer a convenient means of realizing such 
a procedure through a technique analogous to that of a Ram- 
sey separated field experiment [|l6|]. After selecting the de- 
sired region of phase space and preparing atoms in the F = 4, 
rap = state, a 50:50 superposition of the two hyperfine lev- 
els can be created by employing a -| pulse of microwave radi- 
ation resonant with the 9.19 GHz F = 3, m F = — ► F = 4, 
mp = magnetic dipole transition. The QDKR light pulses 
can then be applied using the same laser beam and pulse en- 
velope as the accelerator mode, but with the key change of 
moving the standing wave to make a = [H]. Importantly, 
each of the states in the superposition would experience the 
standing wave potential with a different strength. For exam- 
ple, light detuned 30 GHz below resonance for atoms in the 
F = 4 level is detuned approximately 40 GHz for F = 3 
atoms. Since the light shift scales inversely with detuning, 
the component of the superposition in the F = 3 state sees a 
kicking potential which is ~ 0.75 the size experienced by the 



component in F = 4. The final step of our proposal would 
be to expose the superposition to a second J microwave pulse 
which differs in phase by an amount Acb from the first pulse. 
If the overlap between the two final states is good then a sinu- 
soidal variation in the population of either state would be ob- 
served as Acb were systematically changed. The exact value 
of the overlap could be determined by measuring the popula- 
tion at Act> = and Acb — ? and then using a method similar 
to that outlined in Ref. [[[(J. 

To summarize, we have demonstrated a new method of 
preparing atoms in a narrow region of phase space from an 
initially broad distribution. This was achieved using quantum 
accelerator modes, a process in which atoms are exposed to 
pulses of an off-resonant standing light wave which is accel- 
erating relative to the atoms' frame of reference. Since the 
experimental configuration is very similar to that used for ex- 
periments which have studied the QDKR, it should be rel- 
atively straightforward to apply this technique to studies of 
quantum chaos. Previous experiments have concentrated on 
the ensemble behavior of the atoms. The work presented here 
paves the way for a more detailed examination of the dynam- 
ics. We have also suggested an experiment which would allow 
the production of a quantum phase space stability map. By 
using a method based on Ramsey separated fields we can find 
the amount of overlap between identical initial states which 
have interacted with chaotic potentials of slightly different 
strengths. We hope that the combination of these two ideas 
will eventually lead to a fuller understanding of the concepts 
of quantum chaos. 
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